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Abstract 

We study the effect of site occupation on the description of quan- 
tum spin systems at finite temperature and mean-field level. We im- 
pose each lattice site to be occupied by a single electron. This is 
realized by means of a specific prescription. The outcome of the pre- 
scription is compared to the result obtained by means of a projection 
procedure which fixes the site occupation to one particle per site on 
an average. The comparison is performed for different representa- 
tions of the Hamiltonian in Fock space leading to different types of 
mean- field solutions. The behaviour of order parameters is analyzed 
for each choice of the mean-field and constraint which fixes the occu- 
pation rate at each site. Sizable quantitative differences between the 
outcomes obtained with the two different constraints are observed. 
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1 Introduction 



We consider ordered quantum spin systems at finite temperature in which 
each lattice site is occupied by one electron with a given spin. Such a configu- 
ration can be constructed by means of constraints imposed through a specific 
projection operation pQ which fixes the occupation in a strict sense. The 
constraint can also be implemented on the average by means of a Lagrange 
multiplier procedure [2]. It is the aim of the present work to confront these 
two approaches in the framework of Heisenberg-type models. 

The description of strongly interacting quantum spin systems at finite 
temperature generally goes through a saddle point procedure which is a ze- 
roth order approximation of the partition function. The so generated mean- 
field solution is aimed to provide a qualitatively realistic approximation of 
the exact solution. 

However mean-field solutions are not unique. The implementation of a 
mean-field structure is for a large part subject to an educated guess which 
should rest on essential properties of the considered system, in particular 
its symmetries. This generates a major difficulty A considerable amount 
of work on this point has been made and a huge litterature on the subject 
is available. In particular, systems which are described by Heisenberg-type 
models without frustration are seemingly well described by ferromagnetic or 
antiferromagnetic (AF) Neel states at temperature T = [31 Hj. It may how- 
ever no longer be the case for many systems which are of low dimensionality 
(d < 2) and (or) frustrated jHHUEj- These systems show specific features. 
An extensive analysis and discussion in space dimension d = 2 has recently 
been presented by Wen ^Uj. The competition between AF and chiral spin 
state order has been the object of very recent investigations in the framework 
of continuum quantum field approaches at T = temperature, see [HI EI - 

The reason for the specific behaviour of low dimensional systems may 
qualitatively be related to the fact that low dimensionality induces strong 
quantum fluctuations, hence disorder which destroys the AF order. This 
motivates a transcription of the Hamiltonian in terms of composite opera- 
tors which we call "diffusons" and "cooperons" below, with the hope that 
the actual symmetries different from those which are generated by AF order 
are better taken into account at the mean-field level j2]. 
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In the present work we aim to work out a rigorous vs. average treatment 
of the half-filling occupation constraint on systems governed by Heisenberg- 
type Hamiltonians at the mean-field level, for different types of order. 

The outline of the paper is the following. In section 2 we first recall the 
projection procedure which leads to a rigorous site occupation. Section 3 is 
devoted to the confrontation of the magnetization obtained through this pro- 
cedure with the result obtained by means of an average projection procedure 
in the framework of the mean-field approach characterized by a Neel state. 
The same confrontation is done in section 4 for the order parameter which 
characterizes the system when its Hamiltonian is written in terms of so called 
Abrikosov fermions [HI 12] in d = 2 space dimensions. In section 5 we show 
that the rigorous projection P is no longer applicable when the Hamiltonian 
is written in terms of composite "cooperon" operators. Conclusions, further 
investigations and extensions are presented in section 6. 

2 Site occupation constraint for quantum spin 
systems at finite temperature T. 

Heisenberg quantum spin Hamiltonians of the type 



i,3 

with {Jij} > can be projected onto Fock space by means of the transfor- 
mation 




(1) 
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where {a^, a^} are anticommuting fermion operators. 

This transformation is not bijective because the dimensionality of Fock 
space is larger than the dimensionality of the space in which the spin oper- 
ators {Si} are acting. Indeed, in Fock space, each site i can be occupied by 
0, 1 or 2 fermions corresponding to the states |0, >, |1, >, |0, 1 >, |1, 1 > 
where 1 , > is the particle vacuum, |1,0 >= | + 1/2 > ,|0, 1 >= | — 1/2 > 
and |1,1 >= | + 1/2,-1/2 > in terms of spin 1/2 projections. Since one 
wants to keep states with one fermion per site only states |0, > and |1, 1 > 
have to be eliminated. This can be performed on the partition function for 
a system at inverse temperature (3 

Z = Tre-? H 

where the trace is taken over the whole Hilbert space by adding a projection 
term 

Z = Tre-^ H -» N) (2) 

where N is the particle number operator and fi = in/2/3 an imaginary chem- 
ical potential [Tj. Indeed, the presence of the states |0, >j and |1, 1 >j on 
site i leads in Z to phase contributions which eliminate each other 

e i*o + j** = 

and hence the contributions of these spurious states are cancelled as a whole. 

The common alternative approximate projection procedure would be to 
introduce a chemical potential in terms of real Lagrange multipliers {Aj} and 
to fix this quantity by means of a saddle point procedure with respect to the 
{A,}s 

Z = Tre-^H ! dXie^-^ (4) 
i J 

where is the particle number operator on site i. 

In the following we compare the outcome of these two projection pro- 
cedures in the framework of different mean-field approximations of systems 
described by Heisenberg Hamiltonians. 
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3 Antiferromagnetic mean-field ansatz 
3.1 Exact occupation procedure 

Starting with the Hamiltonian defined in Eq. the partition function Z 
can be written in the form 

z= /n p ({^a.})e- A({ ^' U) (5) 

J i,a 

where the £i 1<T } are Grassmann variables corresponding to the operators 
{a irT , die} defined above. They depend on the imaginary time r in the interval 
[0, 0\. The action A is given by 

A(Ui a ,^}) = f P drJ2^Ar)d T ^ a (r) + W({^(r),^(r)})) (6) 
where 

H{t) = H(r) - hN{t) (7) 

and N(t) is the particle number operator. A Hubbard - Stratonovich (HS) 
transformation which generates the vector fields leads to the partition 
function Z which can be written in the form 

Z = [ II V (i^ €w ^» e ~ ! ° dT[E ''° ^ {t)9 ^ At)+A{{ ^^' ' })] (8) 

In Eq. (jHJ) the expression of Z is quadratic in the Grassmann variables 
{Co-' £«,o-} over which the expression can be integrated. The remaining ex- 
pression depends on the fields {<fi(r)}. A saddle point procedure decomposes 
0i(r) into a mean-field contribution and a fluctuating term 

0,(r)= 0^ + 50^) (9) 
where <^ are the constant solutions of the self-consistent equation 

(10> 
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The partition function takes the form 



Z = Z mS J D({^})e- A «^ (11) 

where the first term on the r.h.s. corresponds to the mean-field contribution 
and the second term describes the contributions of the quantum fluctuations. 

Considering the mean-field contribution the partition function Z m f can 
be put in the form 

Z mf = Tre-^f 

which reads 

Z m/ = r V S 5' J >■"" " n( l + e-^)(l + e -«-) (12) 

where M is the number of sites. The energies Ef and E~ are obtained 
through a diagonalization of the Hamiltonian by means of a Bogolioubov 
transformation in Fock space. Explicitly 



E? = V + " 



W) i 

! ..... 



2 

Working out the expression of Z m f leads to 



and the free energy is given by the expression 

>w = - jDj")^^ - 1 E f^p 1 ) (14) 

ij i \ 
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In order to obtain the expression of the local magnetization {mj one 
should add the term ^ BiSi to the Hamiltonian H in Eq. (0). Going through 
the same steps as above the local magnetization given by 



m 



(mf) 



mf 



8B. 



Bi=0 



(15) 



is related to the {(^"^j's by 



ITLi 



which leads to the self-consistent equation for the {m| m }'s 



rn 



(mf) 



m 



(WO 



(16) 



If the local fields {-Bj} are oriented along a fixed direction e z , m 



(mf) 



(mf) 



e z , the magnetizations are the solutions of the self-consistent equations 



rn 



(mf) 



-th 
2 



(mf)' 



(17) 



3.2 Lagrange multiplier approximation 

Similarly to the preceding case one can introduce the one-particle site oc- 
cupation by means of a Lagrange procedure added on the expression of the 
Hamiltonian 



where A is a variational parameter and {rij} are particle number operators. 
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Following the same lines as above with the help of a HS transformation 
and staying in the ordinary space representation the mean-field partition 
function Z^j can be worked out and reads 

zit = e ^ ; « n( i + e -^) (1+e -^) (is) 



with 



W)i 



W 3 



W) | 



The parameter A is fixed through a minimization of the corresponding 
free energy T^) with respect to this multiplier. The minimization shows 
that the extremum solution is obtained for A = and 

= -\E( J -%^M a " - |E toM (^f^j da) 

which is different from the expression of Eq. (fTlj) . 

The magnetization can be obtained in the same way as done above. One 
obtains 



™aa = \t» { 0Z ' J f^ i (20) 



which is again different from the expression obtained in the case of a rigorous 
projection, see Eq. (JT7J). 

The uniform solutions m\ m ^ = {—l) % rrS m ^ and ra^f^ = (— l)Vn,^o' for 
{Jij} = J have been calculated by solving the selfconsistent equations (17) 
and (20). 

The results are shown in Fig. (fl|) . It is seen that the treatment of the site 
occupation affects sizably the quantitive behaviour of observables. In partic- 
ular it shifts the location of the critical temperature T c by a factor 2. Such a 
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Figure 1: Magnetization vs. reduced temperature t = T/\J\. Full line: exact 
site occupation. Dashed line: average site occupation. 

strong effect has already been observed on the behaviour of the specific heat, 
see refs. [TBI fTOj . 



4 Spin state mean-field ansatz in 2d 

In 2d space the Heisenberg Hamiltonian given by Eq.(l) can be written in 
terms of composite non-local operators {V^} ("diffusons") j2] defined as 

Vij = a^dji + a\ l a jl 
Then the initial Hamiltonian Eq. with J^- = J takes the form 

^ = -^E(H 2? «-y + 2 4 1) (21) 

<ij> 

where i and j are nearest neighbour sites. The physical justification for such 
a transcription relies on the fact that systems with low space dimensionality 
at finite temperature show physical symmetry properties which are closer to 
a spin liquid than an ordered Neel state as already mentioned in the intro- 
duction. 
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The number operator products {riirij} in Eq. (|21|) are quartic in terms of 
creation and annihilation operators in Fock space. In principle its existence 
necessitates the introduction of a further mean field. One can however show 
that its elimination has no influence on the results obtained from the parti- 
tion function. As a consequence we leave it out from the beginning as well 
as the contribution corresponding to the {n^} term. 



4.1 Exact occupation procedure 

Starting with the Hamiltonian 



H = ~J2 V h V v-V N ( 22 ) 

<ij> 

one performs an HS transformation on the corresponding functional inte- 
gral partition function in which the action contains the occupation number 
operator as written out in Eq. (|7J). The Hamiltonian takes the form 

H =TJ\H A 'A.; + E [ A v V v + A v V l] ~ /' V (23) 

' ' <ij> <ij> 

The fields A^ and their conjugates A^ can be decomposed into a mean-field 
contribution and a fluctuation term 

Ay = AT/ + 5A tJ (24) 



The field A™/ can be chosen as a complex quantity A™ - ^ = |A^|e*^ . 
Consider a square plaquette (i, i + e x , i + e y , i + e x + e y ) where e x and e y are 
the unit vectors along the directions Ox and Oy starting from site i on the 
lattice. On this plaquette 

0= n c f 

(ij)epl 

is taken to be constant. If the gauge phase (p™/ fluctuates in such a way 
that 4> keeps constant the average of A™/ will be equal to zero in agreement 
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with Elitzur's theorem ^7j. In order to guarantee SU(2) invariance of the 
mean-field Hamiltonian along the Wilson loop on the plaquette we follow 
HQ QH EE2 EES El and introduce 



e i -I( !) 1 jf f = ? . _)_ g 

^sw. ' " (25) 

e 4 v ^ , it rj = ri + e y 



Then the total flux through the plaquette 



(p — 7T 

In the mean-field approximation the corresponding partition function 
reads 

Z mf = e -P{H mf ^N) ( 26 ) 

where 

H «* = ]§| E A ^- A ^ + E [ A S /p ^ + A ^4] - ^ (27) 

<«i> <«i> 



Performing a Bogolioubov transformation which diagonalizes the remain- 
ing expression in Fourier space leads to 

fc,cr 

where z = 4 is the coordination and A/" the number of sites. The momenta 
{k} act in the first half Brillouin zone (spin Brillouin zone) and the operators 
/?} are the transformed of the {a\ a} 's through the rotation which leads 



to the diagonal expression Eq. (1281) . The eigenenergies Et and E- are 
given by 

= -fi + 2A[cos 2 (k x ) + cos 2 (k y )] 1/2 (29) 

and similarly 

EZ^ = -n- 2A[cos 2 (k x ) + cos 2 (k y )] 1/2 (30) 
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The partition function Z m f has the same structure as the corresponding 
partition function in Eq. (JT5J) and the free energy is given by 

*w = ^ — ^E ln2 ( c ^ Ae *) ( 31 ) 

with 

eg = 2[cos 2 (^) + cos 2 ^)] 1 / 2 (32) 



Finally the variation of T m $ with respect to A leads to the self-consistent 
mean-field equation 

i4E^W (33) 



2JV 

fc,(T 



with A = zA/|J|. 



4.2 Lagrange multiplier approximation 

Similarly to Eq. ([230 one may introduce a Lagrange constraint and write 



H W = T E A <A, + E ( A '/ P - + AyPj) + E " 1) (34) 

<ij> <ij> i 

Then for Aj = A 
with the eigenenergies 

= A + 2A[cos 2 (k x ) + cos 2 (k y )] l/2 (36) 

and similarly 

EZ* = A - 2A[cos 2 (/y + cos 2 (fcj,)] 1/2 (37) 
12 



The expression of the free energy is now given by 

fSj = -JVA + - * £ ln[l + e-" E K][l + e-» E C\ (38) 

k,a 



The minimization of this expression in terms of A delivers the solution 
A = and 

^_^ A + ^_lp ln ( 2cft ^) (39) 



The variation of with respect to A leads to the self-consistent mean- 
field equation equation 

k,a \ / 

with &W = zA/|J|. 

Expressions in Eq. ()39j) and Eq. (}4T)|) should be compared to the expressions 
obtained in Eq. (J31|) and Eq. (J33)) . Fig. 2 shows the behaviour of A for the 
two different treatments of site occupation on the lattice. 



4.3 Magnetization 

We calculate the magnetization in the case of a rigorous one-particle site 
occupation. Starting with the Hamiltonian H given by Eq. (|22|) and adding 
the staggered interaction ^^(-1)'^ of the spins with an homogeneous 
external field h one obtains the mean field expression 

k,cr 
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Figure 2: A vs. reduced temperature t = zT/\J\. Full line: exact site 
occupation. Dashed line: average site ocupation. 



where 

and similarly 
with 



h,k,a r h,k,a 



e: 



h,k,ir 



~V - E h,k,a 



(42) 
(43) 



. h a 2 2 

E h,k,a = [ -J + Ae k 



1/2 



The free energy can now be written 

J\fzA 2 .1 



(44) 



and the magnetization 



m (mf) = 



OF, 



mf 



reads 



m (mf) = 



OA 
~dh 



dh 
2MzA 

h=0 \ \J\ 
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h=Q 



fe,CT 



(45) 



— * 

The sum over k runs over half the Brillouin zone. Going back to Eq. (J33|) 
one sees that mS m f> is equal to zero. 

The same result holds in the case of the average Lagrange multiplier proce- 
dure. 



5 Cooperon mean-field ansatz 

Starting from the Hamiltonian 

H = -\J\ (46) 

<i,j> 

one can introduce a further set of non-local composite operators ("cooper- 
ons") 

dj = a^ctji - o^Ojt (47) 
This leads to the expression 

I^E (4C« + =31) (48) 

<i,j> 

where = J2<j n ia- 



5.1 Exact occupation procedure 

As in the preceding cases it is possible to implement a HS procedure on 
the {C^} and {n^} in such a way that the expresssion of the corresponding 
partition function gets quadratic in the fields {a^,a«j}. The corresponding 
HS fields are {Ty}, and 

h - = A Yl + E + T *A + Tj\ ViVj ) + ^ ^ (49) 

<i,j> <M> ' i 

Introducing the homogeneous fields {T = T^} ,{v = v{\ and integrating over 
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the cooperon fields {Cy} one gets in Fourier space 

Mz\Y\ 2 2M\u\ 2 sr . . t . . 

H = -JJ^ + ^- + Y^-v) a L a k.« (50) 



fe,<T 

az 



+ — luiFaz a ? — Tat a) - ) 



2 

k,cr 



with 

7g = ^(cosfc^. + cosk y ) 

The third term in this expression is complex since /i = iir /2{3. In this rep- 
resentation it is not possible to find a unitary Bogolioubov transformation 
which diagonalizes H . Hence a rigorous implementation of the constraint on 
the particle number per site is not possible. 



Reasons for this situation may be the fact that the Hamiltonian contains 
terms with two particles with opposite spin created or annihilated on the 
same site which is incompatible with the fact that such configurations are 
not allowed in the present scheme. Terms of this type are typical in mean- 
field pairing Hamiltonians which lead to a non-conservation of the number 
of particles of the system. 



5.2 Lagrange multiplier approximation 

If the sites are occupied by one electron in the average the Lagrange proce- 
dure works. We do not develop the derivation of the mean-field here since it 
has been done elsewhere j2]. 



6 Conclusions. 

In summary we have shown that a strict constraint on the site occupation of 
a lattice quantum spin system described by Heisenberg-type models shows 
a sizable quantitative different localization of the critical temperature when 
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compared with the outcome of an average occupation constraint. Conse- 
quently it generates sizable effects on the behaviour of order parameters and 
other physical observables. This is true both in the case of antiferromagnetic 
(Neel) and spin state symmetry. 

Due to the complexity of quantum spin systems the choice of a physically 
meaningful mean field may depend on the coupling strengths of the model 
which describes the systems jTH]. A specific mean-field solution may even 
be a naive way to fix the "classical" contribution to the partition function 
which may in fact contain a mixture of different types of states. As already 
mentioned many efforts have been and are done in order to analyze and over- 
come these problems, see f.i. jHl El ED]- It is our aim to repeat the analysis 
we did above in the framework of a description which is able to implement 
more appropriate descriptions at the mean-field level. 

In a more realistic analysis one should of course take care of the contri- 
butions of quantum fluctuations which may be of overwhelming importance 
particularly in the vicinity of critical points. We expect to work out the first 
order contributions in a loop expansion for both types of constraints and sym- 
metries considered in the present work in order to analyze and compare their 
relative importance as a function of temperature. Such an analysis may pro- 
vide some view on the relative contribution of quantum corrections to order 
parameters and other observables for the different chosen mean- field ansatze. 

The authors would like to thank D. Cabra and T. Vekua for interesting 
and helpful discussions during the time of preparation of the present work. 
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